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Abstract 

The process e+e~ —>■ Z7, where the Z boson decays into hadrons or neutrinos, is 
studied with data collected with the L3 detector at LEP at centre-of-mass energies 
from 189 GeV up to 209 GeV. The cross sections are measured and found to be 
in agreement with the Standard Model predictions. Limits on triple neutral-gauge- 
boson couplings, forbidden in the Standard Model at tree level, are derived. Limits 
on the energy scales at which the anomalous couplings could be manifest are set. 
They range from 0.3 TeV to 2.3 TeV depending on the new physics effect under 
consideration. 
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1 Introduction 



The process e+e~ — > Z7 allows to test the existence of new physics [1], such as anomalous 
couplings between neutral gauge bosons. Effects coming from ZZ7 and Z77 couplings are 
expected to be very small in the Standard Model [1,2], but can be enhanced in compositeness 
models [3] or if new particles enter in higher order corrections. Anomalous ZZ7 and Z77 
couplings would increase the e+e~ Z7 cross section and produce an enhancement of large 
polar angle photons. 

Assuming only Lorentz and U{l)em gauge invariance, the most general form of the ZV^ 
vertices, with a real Z and 7 in the final state, is parametrized by means of the anomalous 
couplings, hf (i = 1 . . . 4; 1/ = 7, Z) [4]. The couplings hY and are CP violating whereas 
and h-Y are CP conserving. All these couplings are zero at tree level in the Standard Model, 
and only the CP conserving ones are non- vanishing (~ 10~"^) at the one-loop level [1,2]. An 
alternative parametrization, which introduces the energy scales of new physics, Ajy, is [5]: 
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where a is the fine-structure constant. The fact that there are always two identical particles at 
the vertex forbids the three bosons to be on-shell. This means the ZZ7 and Z77 vertices may 
only appear if one of the bosons is off-shell. A treatment of these vertices where all three bosons 
are off-shell is discussed in Reference 6. In this Letter, the e+c^ Z7 process is analyzed. 
The maximal experimental sensitivity is achieved with the analysis of the c+c^ ii^ process, 
with the fermion pair in the vicinity of the Z resonance, where the signal statistics is high 
and backgrounds are reduced. In this scenario, effects from an off-shell final-state Z boson are 
negligible [7]. 

L3 published results [8] on the e^e^ Z7 process, setting limits on ZZ7 and Z77 couplings 
from data obtained at lower center-of-mass energies (161 GeV < ^/s < 189 GeV). Results have 
also been pubhshed by other experiments at LEP [9] and at the TEVATRON [10]. In this 
Letter we present results for the highest energies collected at LEP. 

The phase space definition for the e"'"e~ Z7 process requires a photon with energy greater 
than 20 GeV and polar angle in the range 5° < < 175°. Every cross section and acceptance 
in this Letter is referred to this fiducial region . 



2 Data and Monte Carlo Samples 

Data collected by the L3 detector [11] at a/s — 189 GeV — 209 GeV with a total luminosity 
of about 626pb~^ are used to study the e"^e~ — > Z7 process in the channels e"^e~ — > qq7 and 
e"^e~ z/P7. 

The Standard Model processes giving rise to these final states are modelled with KK2f, for 
e'^e" —>■ qq7(7), and KKMC, for e'^'e" — > 1/1/7(7) [12]. Both programs are general purpose 
Monte Carlo generators for the process e'^e" — > ff -1-717, containing complete 0{a'^) corrections 
from initial- and final-state radiation, including their interference. Data for e+e~ — > 1/1/7(7) at 
^/s = 189 GeV, previously analysed and compared to the KORALZ Monte Carlo [13], are now 
re-analysed using KKMC. 
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Background processes are simulated with EXCALIBUR [14] for the four-fermion final states, 
PHOJET [15] and DIAG36 [16] for two-photon collisions with hadrons or Icptons in the final 
state, respectively, and BHWIDE [17] and TEEGG [18] for e+e" e+e--f{-f) . 

All generated events are passed through a simulation of the L3 detector [19] and the same 
analysis procedure as used for the data. Time-dependent detector inefficiencies, monitored 
during data-taking, are also taken into account. 

3 Event Selection 

3.1 Photon Selection 

The main signature of the process e'^'e" — > Z7 is the production of a high energy photon. A 
photon candidate is identified as a shower in the barrel or endcap region of the BGO crys- 
tal electromagnetic calorimeter, consistent with an electromagnetic shower and with a mini- 
mum energy of 5 GeV. The mass of the system recoiling against the photon of energy E^: 
nT'rec — {s — 2£^^y^)^/^, is required to satisfy 80 GeV < m^ec < 110 GeV, consistent with Z- 
boson production. For the y/s values considered, the cuts on the recoiling mass correspond to 
photon energies between 62 GeV and 89 GeV. 

3.2 Selection of e+e" qq7 Events 

In addition to requiring a photon candidate, e"'"e~ — > qq7 events are selected by demanding 
that: 

• the event have more than 6 charged tracks reconstructed in the fiducial volume of the 
tracking chamber and more than 11 calorimetric clusters in the electromagnetic calorime- 
ter. 

• the transverse energy imbalance be less than 15% of the total reconstructed energy and 
the longitudinal energy imbalance less than 20% of the same quantity. 

• the angle of the photon candidate with respect to the beam direction, 9j, satisfy ] cos < 
0.97. 

In order to reject electrons produced in the central region, photon candidates with | cos 9^ \ < 0.90 
are not considered if they are associated to a charged track in the central tracking chamber. 
This requirement eliminates a substantial contamination of background processes. For the last 
period of data taking, corresponding to data at > 202 GeV, this rejection cut is relaxed to 
I cos 6'^ I < 0.85 to account for different running conditions of the detector. This change increases 
the contamination from background processes. 

Table 1 hsts the data luminosity analysed at each ^/s, the selection efficiency, the background 
level and the number of selected events, after background subtraction. 

The trigger inefficiency is estimated to be negligible due to the redundancy of subtriggers 
involved in tagging this final state. Two backgrounds contribute in equal proportions: the 
e'^e" — > qq'ev and e"^e~ — > qqe'^e" processes, where one electron fakes a photon. 

Figure 1 shows the distributions of TTZyg^^ COS 9 J and the invariant mass of the hadron system, 
reconstructed from jet and photon directions and ^/s. 

The resolution of the L3 electromagnetic calorimeter, better than 1%, allows the observation 
of a tail in Figure la, for rrirec above the nominal Z mass, due to initial state radiation photons. 
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3.3 Selection of e+e i/i/j Events 

In addition to the presence of a photon, selected as described above, the events from the 
e+e~ — > i/P^ process are selected by the following criteria: 

• the event must have at most 5 calorimetric clusters, due to low energy (< 1 GeV) initial 
state photons or noise in the calorimeter. The number of hits in the tracking chamber 
associated to a calorimetric cluster, including the photon candidate, must not exceed 40% 
of the expected number of hits for a charged track. 

• the angle of the photon candidate with respect to the beam direction must satisfy | cos 9j\ < 
0.96. This range differs from that of the hadronic channel in order to match the angular 
coverage of the central tracking chamber used to reject electrons. 

• the total reconstructed energy, Etot, must fulfill ^/s — Efot > 0.95 Etot and the transverse 
energy imbalance must be greater than 0.2Etot- 

• To suppress cosmic ray background, there must be at least one scintillator time measure- 
ment within ±5 ns of the beam crossing time. The scintillator signals must be associated 
with calorimetric clusters. 

The background in the selected sample is found to be negligible. Table 2 lists the data lu- 
minosity analysed at each a/s and the selection efficiency. The selection efficiency includes 
trigger efficiency, evaluated to be around 95% by using two independent data samples from the 
e"'"e~ — > e"'"e~ and c+c~ — > 77 processes. The number of selected events is also given. Figure 2 
shows the distributions of m^ec and cos 6'^. 

4 Cross Section Measurements 

The measured cross sections for both the e+e~ qqj and e+e" — > uP^ processes are presented 
in Tables 1 and 2, respectively, together with Standard Model predictions. Good agreement 
is observed. The uncertainty on the expected cross section, asM, takes into account a 1% 
theory uncertainty of KK2f and KKMC and the finite Monte Carlo statistics generated for 
these studies. 

In addition to cross sections, in Table 3 we present more detailed information for the e''"e~ — > 
qq7 process on the number of events observed and expected, the background and selection 
efficiencies, in bins of rrirec and | cos^-y|. Similar tables are given elsewhere [20] for the e+e" — > 
uu'-f process. 

The main sources of systematic uncertainties are summarized in Table 4. The largest con- 
tribution is due to the selection procedure. A change of 3% in the values of the cut on rrirec 
corresponds to 0.8% and 1.5% uncertainties for hadronic and invisible decay modes, respectively. 
Changes in the photon energy scale give uncertainties of 0.4% and 0.6% for the hadronic and 
invisible channels, respectively. The uncertainty from limited Monte Carlo statistics amounts 
to 0.4% for the e"'"e~ — >• uuj channel and varies between 0.1% and 0.4% for the e"'"e~ — >• qq7 
channel. The accuracy on the luminosity estimation gives a 0.2% uncertainty. Uncertainty in 
the measurement of the trigger efficiency contributes an additional 0.3% in the e"'"e~ z/z/7 
process. A variation of 10% in the background level corresponds to a 0.3% uncertainty in the 
hadronic channel. 
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The variation of the sum of the e''"e~ — > qq7 and e"'"e~ —>■ uu^ cross sections with y/s 
is presented in Figure 3. Cross sections at = 161, 172 and 183 GeV, aheady pubhshed 
by L3 [8], are included for completeness. The relative deviation from the Standard Model 
predictions as a function of y/s is also shown. Good agreement is found. 



5 Triple Neutral-Gauge-Boson Couplings 

Since deviations from Standard Model expectations are found neither for the e'*'e~ — > qq7 
nor for the e''"e~ — > i/Pj process, limits on anomalous triple-neutral-gauge boson couplings are 
extracted by using an optimal observable method [21]. 



5.1 Optimal Observable Method 

In the presence of anomalous couplings, the cross section for the process e"'"e~ — > Z7 is propor- 
tional to \MsM + MacI"^, with MsM and MAcihJ) the Standard Model and anomalous coupling 
amplitudes, respectively. 

As Mac depends hnearly on the hf (i=l . . .4; V=7,Z) parameters, the differential cross 
section can be written as a quadratic function on the anomalous couplings: 

^ = co{n) + j: E c..Am + i: E E E (3) 

— * — * 

where Q stands for the phase space variables defining the final state, Co(fi) is the Standard 
Model cross section and Ci ^ ^ and Ca.ij.v-.v' are coefficients related to the anomalous amplitudes. 
The variables defined for each coupling as: 

o.,,vm ^ (4) 



02,ij,v,v'{^l) = — (5) 

are called "optimal variables" as they contain the full kinematic information on the event and 
allow the determination of the parameters hf with the maximum possible statistical precision. 
If the parameters hf are sufficiently small, the quadratic term can be neglected and all the 

— * 

information in the multidimensional phase space is projected into the variable Oi. 

In order to extract the hj^s, a binned maximum-likelihood fit of the expected distribution 
of the optimal variable Ci,i,y is performed to the data, assuming a Poisson density distribution 
in each bin. Both the shape of the optimal variable distribution, which includes energy and 
angular information, and the total number of events contribute to the fit. The expected number 
of events in the presence of anomalous couplings is computed from a Standard Model reference 
sample by applying a rcwcighting technique, in which each Monte Carlo event is weighted with 
the following quantity, defined at the generator level: 

\JyisMr 
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The comparison between expected and observed events is done at the level of reconstructed 
variables so that all experimental effects, such as detector resolution or selection efficiencies, 
are automatically taken into account. 

5.2 Limits on Anomalous Couplings 

Making use of the optimal observable method and taking into consideration the information 
on the total event rate for each process and the phase space variables defining the final state, 
limits at 95% Confidence Level (CL) are set on the couplings. The reconstructed set of 
variables used to compute the optimal variables is Q = {E^, 9.y, (p^, 9^, 0^), where E-y, 6^ and 0-y 
are the energy and angles of the photon, and 6'^ and 0^ the angles of the fermion / in the Z 
rest frame. In the e+e" — > vv^ channel only the three photon variables are used. 

Distributions of the optimal variables for the couplings and /i] are shown in Figure 4. The 
regions of maximal sensitivity to the existence of anomalous couplings correspond to the largest 
absolute values of the optimal variables, where discrepancies with Standard Model predictions 
are expected to be larger. 

The 95% CL limits on each individual anomalous coupling, combining both channels and 
from all data collected at 189 GeV < < 209 GcV, are obtained from one-dimensional fits. 
The results are given in Table 5 and they correspond to the following intervals: 

-0.153 <hl< 0.141 -0.057 < Kl < 0.057 

-0.087 <hl< 0.079 -0.050 < hi < 0.023 

-0.220 < /if < 0.112 -0.059 < Z^;^ < 0.004 

-0.068 <hl< 0.148 -0.004 < hj < 0.042 

To obtain these intervals one parameter is left free at a time, setting the other seven anoma- 
lous couplings to zero. These limits supersede the previous L3 results, obtained with a smaller 
data sample at lower centre-of-mass energies [8]. The observed limits agree within 10% with 
the expected limits. 

Limits coming from even couplings, and h\ are more stringent than those coming from 
odd couplings, h\ and /i^, as the former correspond to operators of dimension eight, while the 
latter correspond to dimension six operators [1,2], as reflected in equations (1) and (2) in the 
different dependence of the parameters on the mass mz and the energy scales. These relations 
imply that the Standard Model is tested more stringently in the linear expansion of the effective 
Lagrangian when considering the even couplings. 

Fits to the two-dimensional distributions of the optimal observables are performed to de- 
termine the pairs of CP-violating {h\,h\) and CP-conserving couphngs {hX,,h\), keeping in 
each case the other couplings fixed at zero. Results at 95% CL are shown in Table 6. A strong 
correlation between the two CP-violating or CP-conserving parameters is observed. Contours 
for the 68% and 95% CL two-dimensional limits on each pair of couplings are shown in Figure 
5. The main sources of systematic uncertainties, discussed in section 4, are included in the 
limits calculation. They contribute 0.02 to one-dimensional limits and 0.03 for two-dimensional 
limits. 

If the data are interpreted in terms of new physics scales using formulae (1) and (2), lower 
hmits at 95% CL on the scale of new physics are obtained as: 
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Alz > 0.8 TeV Ai^ > 1.3 TeV 

Asz > 0.3 TeV Aa^ > 0.4 TcV 

Asz > 0.8 TeV As^ > 2.3 TeV 

A4Z > 0.3 TeV K^^ > 0.4 TeV. 

To determine the confidence levels the probability distributions are normalized over the 
physical range of the parameters, A > 0. 



6 Conclusions 

The analysis of the process e+e~ — > Z7 in the final states Z7 — > qq7 and Z7 — > 1/P7 with 
620 pb~^ of luminosity collected by the L3 detector at 189 < ^/s < 209 GeV reveals a very good 
agreement between the measured cross sections and the Standard Model prediction. Detailed 
information on the hadronic final-state events are given in form of tables to allow constraints 
of future models. 

These measurements establish upper limits at 95 % CL on the values of anomalous couplings, 
hj , appearing in the triple neutral-gauge-boson vertices, ZZ7 and Z77. At tree level in the 
Standard Model these couplings are zero. We observe no deviation from this prediction and 
constrain possible values of the anomalous couplings in intervals of widths between 0.05 and 
0.33, depending on the coupling considered. These limits improve and supersede our previous 
limits [8]. 
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(GoV) 




r (<X) 


Back. ('X) 


E\xnits 


a (pi)) 


(^SM (pb) 


188.6 


172.1 


28.9 ± 0.1 


0.9 


899 


18.1 ± 0.6 ± 0.2 


18.7 ± 0.2 


191.6 


17.9 


28.2 ± 0.4 


0.8 


101 


20.0 ± 2.1 ± 0.2 


18.8 ± 0.2 


195.5 


74.9 


27.1 ± 0.2 


0.9 


351 


17.3 ± 1.1 ± 0.2 


17.9 ± 0.2 


199.5 


67.4 


27.3 ± 0.2 


1.0 


333 


18.1 ± 1.0 ± 0.2 


16.9 ± 0.2 


201.7 


36.5 


27.4 ± 0.3 


1.0 


163 


16.3 ± 1.3 ± 0.2 


16.3 ± 0.2 


202.5 - 205.5 


78.7 


25.0 ± 0.1 


4.8 


325 


16.4 ± 1.0 ± 0.2 


16.3 ± 0.2 


205.5 - 207.2 


124.1 


25.2 ± 0.1 


4.7 


494 


15.8 ± 0.7 ± 0.2 


16.1 ± 0.2 


207.2 - 209.2 


8.2 


24.8 ± 0.1 


4.7 


33 


16.1 ± 3.0 ± 0.2 


15.6 ± 0.2 



Table 1: Integrated luminosities, L, and results of the e+e" qcp/ selection: se- 
lection efficiencies, e, background level, number of selected events (background sub- 
tracted) and measured cross sections with statistical and systematic uncertainties. 
The uncertainty on e comes from Monte Carlo statistics. The corresponding Stan- 
dard Model cross sections, asu, are listed in the last column. They are derived 
from the KK2f Monte Carlo generator [12]. Their uncertainty includes a 1% theory 
uncertainty and finite Monte Carlo statistics. 



(GeV) 


C (pb-i) 


e(%) 


Events 


a (pb) 


CTSM (pb) 


188.6 


175.6 


32.3 ± 0.4 


288 


5.1 ± 0.3 ± 0.1 


4.99 ± 0.05 


191.6 


16.9 


31.4 ± 0.5 


25 


4.7 ± 0.8 ± 0.1 


4.85 ± 0.05 


195.5 


80.9 


31.5 ± 0.4 


107 


4.2 ± 0.4 ± 0.1 


4.57 ± 0.05 


199.5 


79.5 


28.0 ± 0.4 


98 


4.4 ± 0.4 ± 0.1 


4.42 ± 0.05 


201.7 


36.1 


30.8 ± 0.4 


50 


4.5 ± 0.6 ± 0.1 


4.34 ± 0.04 


202.5 - 205.5 


74.3 


28.9 ± 0.5 


88 


4.1 ± 0.4 ± 0.1 


4.17 ± 0.04 


205.5 - 209.2 


129.6 


29.4 ± 0.5 


160 


4.2 ± 0.3 ± 0.1 


4.09 ± 0.04 



Table 2: Integrated luminosities, £, and results of the e+e" — > selection: selection efficien- 
cies, e, number of selected events and measured cross sections with statistical and systematic 
uncertainties. The uncertainty on e comes from Monte Carlo statistics. The corresponding 
Standard Model cross sections, ctsm, are listed in the last column. They are derived from the 
KKMC Monte Carlo generator [12]. Their uncertainty includes a 1% theory uncertainty and a 
contribution from finite Monte Carlo statistics. 
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COS 9-. 



0.00-0.20 
0.20 - 0.40 
0.40 - 0.60 
0.60 - 0.80 
0.80 - 0.90 
0.90 - 0.99 



80-88 



88-92 



[GeV] 
92-96 



96 - 104 



104 - 110 



14/20.4/0.0/78 
19/27.1/0.0/85 
26/25.7/0.0/62 
29/19.3/0.2/28 
35/39.7/0.6/61 
95/55.9/4.6/18 



64/ 69.8/ 0.0/35 79/ 62.3/ 0.0/35 

91/ 86.4/ 0.0/37 92/ 75.4/ 0.9/38 

116/110.1/ 0.4/37 121/107.9/ 0.5/41 

88/ 84.1/ 0.5/17 75/100.6/ 0.1/24 

145/167.4/ 0.4/35 127/159.6/ 1.9/40 

302/276.8/11.5/12 271/281.9/11.8/15 



58/ 55.8/ 0.5/57 43/24.0/0.3/67 

62/ 68.2/ 0.3/62 31/26.9/0.7/64 

79/ 83.8/ 0.6/64 35/32.9/0.6/64 

82/ 92.3/ 1.3/43 37/40.7/0.5/48 

161/136.1/ 4.3/63 68/55.7/1.9/70 

232/225.8/12.3/24 98/87.1/4.4/24 



Table 3: Numbers of events selected in the e+e~ qq7 channel, Standard Model 
expectations, background level and selection efficiencies (in %) as a function of the 
recoil mass, rrirec to the photon candidate and of the absolute value of the cosine of 
the photon polar angle, | cos 9^ \ . 





Uncertainty (%) 


Source 


e"'"e — s> qq7 


e^e — > 


Selection Process 


0.8 


1.5 


Photon energy scale 


0.4 


0.6 


MC statistics 


< 0.4 


0.4 


Luminosity 


0.2 


0.2 


Trigger efficiency 




0.3 


Background level 


0.3 




Total 


1.1 


1.7 



Table 4: Sources of systematic uncertainty in the e"'"e — > qq7 and e"'"e — > vv^ 
cross sections. 





Fitted 


Negative 


Positive 


Parameter 


value 


error (95% CL) 


error (95% CL) 


hi 


-0.007 


0.146 


0.148 


hi 


-0.006 


0.080 


0.085 


hi 


-0.036 


0.184 


0.148 


hi 


0.038 


0.106 


0.110 


hi 


-0.001 


0.056 


0.058 


hi 


-0.014 


0.035 


0.037 


hi 


-0.026 


0.033 


0.031 


hi 


0.020 


0.024 


0.022 



Table 5: Fitted values and errors at 95 % CL on the individual anomalous couphngs 
from one-dimensional fits. In each fit the other seven parameters are set to zero. 



Parameter 


Fitted 

value 


Negative 
hmits 


Positive 
hmits 


Correlation 
coeficient 


hi 
hi 


-0.05 
-0.03 


-0.38 
-0.22 


0.30 
0.18 


0.89 


hi 
hi 


-0.00 
0.04 


-0.46 
-0.24 


0.40 
0.28 


0.90 


hi 
hi 


-0.04 
-0.03 


-0.15 
-0.09 


0.07 
0.04 


0.85 


hi 
hi 


0.03 
0.04 


-0.09 
-0.04 


0.14 
0.11 


0.93 



Table 6: Limits at 95 % CL on pairs of anomalous couplings from two-dimensional 
fits. In each fit the other six parameters are set to zero. 
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Figure 1: Distributions of a) the recoil mass to the photon candidate in e+e~ qqj 
events b) the polar angle of the photon and c) the invariant mass of the hadron 
system. The points are data, the open histogram is the Standard Model Monte 
Carlo prediction and the hatched one is the background from the e"'"e~ qq'eu and 
e+e~ — > qqe+e~ processes. 
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Figure 2: Distributions of a) the recoil mass to the photon candidate in e+e~ — > vV^ 
events and b) the photon polar angle. The points are data and the histogram is the 
Standard Model Monte Carlo prediction. 
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Figure 3: Variation of the sum of the cross-sections of the e"'"e~ — > qq7 and e"'"e~ — > 
vv"^ processes, 0"(qq7 + 1^2/7), with ^/s. Data are represented by the dots, while the 
sohd hne gives the variation of the Standard Model cross-section, ctsm, as calculated 
with the KK2f and KKMC [12] Monte Carlo programs. The width of the band 
takes into account a 1% uncertainty in each of the e"'"e~ qq7 and e"''e~ uu'f 
theoretical cross sections. The lower plot shows the relative difference. 
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Figure 4: Distribution of the optimal variables for a) the CP-violating coupling and b) for the CP-conserving coupling hj. Data are 
shown together with the expectations for the Standard Model (SM) and for values of anomalous couplings = ±0.5 and hJ = ±0.2. 
The lower plots shows the ratios between the anomalous coupling contributions and the data, to the Standard Model expectation. 
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Figure 5: Two dimensional limits at 68% and 95% CL on the pairs of CP-violating 
coupling parameters, a) /if vs. and b) hi vs. hj and the pairs of CP-conserving 
coupling parameters, c) /if vs. /if and d) hJ vs. hJ. The Standard Model predictions 
are indicated by the points. The shaded areas correspond to the regions allowed at 
68% CL while the dashed line shows the 95% CL exclusion contour. 
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